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Soller  BR,  Yang  Y,  Soyemi  OO,  Ryan  KL,  Rickards  CA,  Walz 
JM,  Heard  SO,  Convertino  VA.  Noninvasively  determined  muscle 
oxygen  saturation  is  an  early  indicator  of  central  hypovolemia  in 
humans.  J  Appl  Physiol  104:  475-481,  2008.  First  published  Novem¬ 
ber  15,  2007;  doi:  10.1 152/japplphysiol. 00600.2007. — Ten  healthy 
human  volunteers  were  subjected  to  progressive  lower  body  negative 
pressure  (LBNP)  to  the  onset  of  cardiovascular  collapse  to  compare 
the  response  of  noninvasively  determined  skin  and  fat  corrected  deep 
muscle  oxygen  saturation  (SmCL)  and  pH  to  standard  hemodynamic 
parameters  for  early  detection  of  imminent  hemodynamic  instability. 
Muscle  SmOj  and  pH  were  determined  with  a  novel  near  infrared 
spectroscopic  (NIRS)  technique.  Heart  rate  (HR)  was  measured  con¬ 
tinuously  via  ECG,  and  arterial  blood  pressure  (BP)  and  stroke 
volume  (SV)  were  obtained  noninvasively  via  Finometer  and  imped¬ 
ance  cardiography  on  a  beat-to-beat  basis.  SmC>2  and  SV  were 
significantly  decreased  during  the  first  LBNP  level  (—15  mmHg), 
whereas  HR  and  BP  were  late  indicators  of  impending  cardiovascular 
collapse.  SmCL  declined  in  parallel  with  SV  and  inversely  with  total 
peripheral  resistance,  suggesting,  in  this  model,  that  SmCL  is  an  early 
indicator  of  a  reduction  in  oxygen  delivery  through  vasoconstriction. 
Muscle  pH  decreased  later,  suggesting  an  imbalance  between  delivery 
and  demand.  Spectroscopic  determination  of  SmCL  is  noninvasive  and 
continuous,  providing  an  early  indication  of  impending  cardiovascular 
collapse  resulting  from  progressive  reduction  in  central  blood  volume. 

tissue  oxygen;  near  infrared  spectroscopy;  physiological  monitoring; 
hemodynamic  instability;  lower  body  negative  pressure 


trauma  and  severe  hemorrhage  remain  the  most  frequent 
causes  of  death  in  the  1-44  yr  age  group  in  both  the  civilian  and 
military  settings  (2,  4,  24).  Early  recognition  of  the  severity  of 
shock  and  prompt  institution  of  appropriate  resuscitative  measures 
are  widely  believed  to  improve  outcome  and  decrease  the  pro¬ 
gression  to  multi-system  organ  failure  by  virtue  of  maintaining 
end-organ  perfusion.  Routine  clinical  parameters  such  as  arterial 
blood  pressure,  heart  rate  (HR),  mental  status,  and  urine  output,  as 
well  as  hematocrit  levels,  are  often  late  indicators  of  the  true 
extent  of  metabolic  derangement  during  hypovolemic  shock.  In  a 
study  on  HR  variability  and  its  association  with  mortality  in 
prehospital  trauma  patients,  Cooke  and  Convertino  (8)  demon¬ 
strated  that  HR,  arterial  oxygen  saturation  (via  pulse  oximetry), 
and  arterial  blood  pressure  were  poor  predictors  of  outcome. 
Hence,  there  exists  a  need  for  noninvasive,  continuous  monitoring 
of  parameters  that  can  provide  an  early  indication  of  central 
hypovolemia. 

The  physiological  response  to  hemorrhage  includes  signifi¬ 
cant  vasoconstriction  to  help  maintain  adequate  perfusion  pres¬ 
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sure  to  vital  organs.  In  swine  models  of  hemorrhage,  this  can 
be  observed  as  an  early,  rapid,  and  significant  decrease  in 
peripheral  and  splanchnic  oxygenation  (3,  26)  as  blood  is 
shunted  to  the  heart  and  brain.  On  the  basis  of  this  known 
physiological  response,  we  hypothesize  that  a  noninvasive 
assessment  of  muscle  oxygen  could  provide  a  significantly 
earlier  indication  of  blood  volume  loss  than  the  standard  vital 
sign  measurements. 

Near  infrared  spectroscopy  (NIRS)  is  in  wide  clinical  use  for 
the  determination  of  arterial  oxygen  saturation  (Sao,),  which 
provides  information  on  pulmonary  gas  exchange.  When  blood 
flow  to  the  peripheral  muscles  is  decreased,  an  increase  in 
oxygen  extraction  is  reflected  in  a  reduction  in  muscle  oxygen 
saturation  (SmCh;  Ref.  3)  and  muscle  oxygen  tension  (Pm02; 
Ref.  28),  both  of  which  can  be  determined  noninvasively  with 
NIRS.  It  is  also  possible  to  determine  muscle  pH  (pHm)  using 
NIRS  (28,  37).  In  swine  hemorrhage,  tissue  pH  is  significantly 
more  sensitive  to  shock  than  measures  of  arterial  and  venous 
pH  (26),  and  depressed  tissue  pH  is  strongly  associated  with 
negative  outcomes  (27).  From  a  practical  standpoint,  NIRS 
instruments  can  be  made  small  and  portable,  allowing  them  to 
be  used  for  rapid,  noninvasive  patient  assessment  both  inside 
and  outside  the  hospital. 

Human  models  for  the  controlled  study  of  acute  hemorrhage 
and  hemorrhagic  shock  are  very  limited  because  of  ethical 
issues  associated  with  the  withdrawal  of  large  volumes  of 
blood  from  healthy  human  subjects  (18).  However,  a  noninva¬ 
sive  human  model  for  assessing  acute,  progressive  central 
hypovolemia  such  as  that  seen  in  hemorrhagic  shock  is  based 
on  the  use  of  lower  body  negative  pressure  (LBNP).  Human 
volunteers  can  be  subjected  to  increasing  levels  of  negative 
pressure  applied  to  the  lower  body,  thus  creating  central 
hypovolemia  that  induces  hemodynamic  responses  similar  to 
those  reported  during  acute  hemorrhage  (9).  We  thus  used  the 
LBNP  model  to  assess  the  response  of  noninvasively  deter¬ 
mined  SmCL  and  pHm  to  central  hypovolemia.  On  the  basis  of 
previous  animal  studies,  we  hypothesized  that  Sm02  and  Pm02 
would  provide  a  very  early  indicator  of  central  hypovolemia, 
while  significant  changes  in  pHm  would  only  occur  later  in  the 
hypovolemic  challenge  when  an  imbalance  between  oxygen 
delivery  and  demand  developed. 

MATERIALS  AND  METHODS 

Subjects.  Ten  healthy,  non-smoking,  normotensive  men  and 
women  (5/5)  with  a  mean  (±SE)  age  of  35  ±  3  yr,  height  173  ±  5  cm, 
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and  weight  75  ±  6  kg  volunteered  to  participate  as  subjects  for  this 
investigation.  All  procedures  and  risks  associated  with  the  study  were 
explained  and  voluntary  written  informed  consent  was  obtained  from 
the  participants.  A  complete  medical  history  and  physical  examination 
was  performed  on  each  of  the  potential  subjects.  In  addition,  female 
subjects  underwent  a  urine  test  to  ensure  that  they  were  not  pregnant. 
Subjects  maintained  their  normal  sleep  pattern,  refrained  from  exer¬ 
cise,  and  abstained  from  caffeine  and  other  autonomic  stimulants  such 
as  prescription  or  non-prescription  drugs  for  at  least  48  hr  prior  to  the 
experimental  protocol  unless  cleared  by  the  physician  medical 
screener  to  continue  taking  the  medications.  All  experimental  proce¬ 
dures  and  protocols  were  reviewed  and  approved  by  the  Institutional 
Review  Board  for  the  use  of  human  subjects  at  the  Brooke  Army 
Medical  Center  at  Fort  Sam  Houston,  TX. 

Protocol.  With  the  use  of  a  neoprene  skirt  designed  to  form  an 
airtight  seal  between  the  subject  and  the  chamber,  the  application  of 
negative  pressure  to  the  lower  body  (below  the  iliac  crest)  results  in  a 
redistribution  of  blood  away  from  the  upper  body  (head  and  heart)  to 
the  lower  extremities  and  abdomen.  This  model  provides  conditions  of 
controlled,  experimentally  induced  hypovolemic  hypotension,  offer¬ 
ing  a  valuable  method  for  investigating  monitoring  devices  such  as 
NIRS.  Although  absolute  equivalence  between  the  magnitudes  of 
negative  pressure  applied  and  actual  blood  loss  cannot  be  determined 
at  this  time,  review  of  human  and  animal  data  has  revealed  ranges  of 
effective  blood  loss  (or  fluid  displacement)  caused  by  LBNP.  On  the 
basis  of  the  magnitude  of  induced  central  hypovolemia,  we  estimate 
that  10-20  mmHg  negative  pressure  produces  hemodynamic  responses 
equivalent  to  those  resulting  from  blood  loss  of  400-550  ml;  20-40 
mmHg  LBNP  induces  hemodynamic  responses  equivalent  to  blood  loss 
of  550-1,000  ml;  while  >40  mmHg  LBNP  induces  responses  equivalent 
to  blood  loss  approximating  1,000  ml  or  more  (9). 

Each  subject  reported  to  the  laboratory  for  a  progressive  LBNP 
protocol  that  was  designed  to  test  his/her  tolerance  to  experimentally 
induced  hypotensive  hypovolemia.  The  subject  was  first  instrumented 
with  noninvasive  devices  to  measure  HR,  stroke  volume  (SV),  arterial 
blood  pressure,  PmCL,  and  pHm.  The  LBNP  protocol  consisted  of  a 
5-min  baseline  period  followed  by  5  min  of  chamber  decompression 
to  — 15,  —30,  —45,  and  —60  mmHg  and  additional  increments  of  — 10 
mmHg  every  5  min  until  either  the  onset  of  cardiovascular  collapse  or 
the  completion  of  5  min  at  — 100  mmHg.  Cardiovascular  collapse  was 
defined  by  one  or  a  combination  of  the  following  criteria:  1 )  a 
precipitous  fall  in  systolic  blood  pressure  (SBP)  >15  mmHg  and/or 
sudden  bradycardia;  2)  progressive  diminution  of  SBP  <70  mmHg;  or 
5)  voluntary  subject  termination  due  to  discomfort  from  presyncopal 
symptoms  such  as  sweating,  nausea,  dizziness,  or  gray-out.  At  the 
onset  of  cardiovascular  collapse,  the  chamber  vacuum  was  released  to 
ambient  pressure  to  rapidly  restore  blood  flow  to  the  central  compart¬ 
ment.  To  assure  subject  safety,  an  ACLS-certified  physician  was 
present  in  the  laboratory  building  during  all  LBNP  tests. 

Hemodynamic  measurements.  Continuous  HR  was  measured  from 
a  standard  electrocardiogram  (ECG).  Beat-by-beat  SBP  and  diastolic 
blood  pressure  (DBP)  were  measured  noninvasively  using  an  infrared 
finger  photoplethysmograph  (Finometer  Blood  Pressure  Monitor, 
TNO-TPD  Biomedical  Instrumentation,  Amsterdam.  The  Nether¬ 
lands).  The  Finometer  blood  pressure  cuff  was  placed  on  the  middle 
finger  of  the  left  hand,  which,  in  turn,  was  laid  at  heart  level.  Excellent 
estimates  of  directly  measured  intra-arterial  pressures  during  various 
physiological  maneuvers  have  been  demonstrated  with  this  device 
(15,  16,  21,  33).  Mean  arterial  pressure  (MAP)  was  calculated  by 
dividing  the  sum  of  SBP  and  twice  DBP  by  three.  Pulse  pressure  (PP) 
was  calculated  by  subtracting  DBP  from  SBP.  Arterial  oxygen  satu¬ 
ration  was  measured  using  pulse  oximetry  (BCI  Capnocheck  Plus; 
Smiths  Medical,  Waukesha,  WI). 

While  the  Finometer  is  capable  of  determining  SV  and  total 
peripheral  resistance  (TPR),  superior  accuracy  in  this  setting  has  been 
observed  with  the  thoracic  bioimpendance  technique  (20).  Beat-to- 
beat  SV  was  measured  noninvasively  using  thoracic  electrical  bioim¬ 


pedance  with  an  HIC-2000  Bio-Electric  Impedance  Cardiograph  (Bio- 
Impedance  Technology,  Chapel  Hill,  NC).  This  technique  is  based  on 
the  resistance  changes  in  the  thorax  to  a  low-intensity  (4  rnA), 
high-frequency  (70  kHz)  alternating  current  applied  by  band  elec¬ 
trodes  placed  at  the  root  of  the  neck  at  the  xiphoid  process  at  the 
midaxillary  line.  Ventricular  SV  was  determined  via  the  Kubicek 
equation  (17):  SV  (in  ml)  =  p  X  (L/Zo)2  X  LVET  X  (dZ/df),  where 
p  (in  ohm/cm)  is  the  blood  resistivity,  a  constant  of  135  ohms/cm  in 
vivo;  L  (in  cm)  is  the  mean  distance  between  the  inner  band  electrodes 
(front  and  back);  Zo  (in  ohms)  is  the  average  thoracic  background 
impedance;  LVET  (in  seconds)  is  the  left  ventricular  ejection  time; 
and  (dZ/df)  is  the  maximum  height  of  the  dZ/d t  peak  measured  from 
the  zero  line.  Correlation  coefficients  of  0.70-0.93  have  been  re¬ 
ported  in  SV  measurements  simultaneously  made  with  thoracic  elec¬ 
trical  bioimpedance  and  thermodilution  techniques  (20).  Cardiac  out¬ 
put  (Q)  was  calculated  as  the  product  of  HR  and  SV,  and  TPR  was 
estimated  by  dividing  MAP  by  Q. 

Hemodynamic  data  were  sampled  at  500  Hz  and  recorded  directly 
to  data-acquisition  software  (WINDAQ,  Dataq  Instruments,  Akron, 
OH).  Analysis  of  data  was  subsequently  accomplished  using  commer¬ 
cially  available  analysis  software  (WinCPRS,  Absolute  Aliens,  Turku, 
Finland).  Data  presented  for  each  of  these  parameters  represent  the 
average  values  taken  over  the  last  3  min  of  baseline  and  each  LBNP 
level. 

Noninvasive  determination  of  SmO 2,  Pm02,  and  pHm.  S1T1O2, 
P1T1O2,  and  pHm  were  determined  noninvasively  using  a  NIRS  mon¬ 
itor  developed  jointly  by  personnel  from  the  Anesthesiology  Depart¬ 
ment  of  the  University  of  Massachusetts  Medical  School  (Worcester, 
MA)  and  Luxtec  (West  Boylston,  MA).  This  spectroscopic  technique 
was  validated  for  individual  heart  surgery  patients  against  invasive 
PO2  and  pH  sensors  inserted  into  the  hypothenar  muscle  (28).  The 
NIRS  system  used  in  this  study  employs  additional  mathematical 
preprocessing  techniques  to  correct  spectra  for  variation  in  skin 
pigmentation,  fat,  and  muscle  optical  properties  prior  to  the  calcula¬ 
tion  of  Sm02,  Pm02,  and  pHm.  These  techniques  are  necessary  to 
allow  one  calibration  equation  for  each  parameter  to  be  used  for  all 
subjects.  The  optical  sensor  collects  NIR  reflectance  spectra  from 
deep  within  the  forearm  muscle  (flexor  digitorum  profundus)  every 
~20  s.  The  spectra  are  then  processed  with  calibration  equations 
contained  in  a  dedicated  computer.  SmCL,  P1T1O2,  and  pHm  are 
simultaneously  calculated  from  each  spectrum,  displayed  as  a  trend, 
and  stored  on  a  hard  drive  contained  in  the  system. 

Light  is  collected  with  two  sensors  contained  in  the  same  housing. 
One  sensor  collects  light  that  illuminates  only  the  skin  and  fat  layer. 
The  second  sensor  collects  light  that  illuminates  the  skin,  fat,  and 
muscle  layer.  Mathematical  processing  removes  the  light  reflected 
from  the  skin  and  fat,  leaving  only  the  absorption  spectrum  of  muscle 
(34).  Removal  of  spectral  interference  from  skin  pigmentation  and  fat 
is  critical  to  determining  absolute  chemical  concentrations  from  mus¬ 
cle  spectra.  The  mathematical  equations  for  determining  S1T1O2  and 
Pm02  from  the  corrected  spectra  are  briefly  described  in  the  appendix. 
A  detailed  description  of  the  SmCL  calculation  and  its  validation  has 
recently  been  published  by  Yang  et  al.  (36).  PmCL  is  calculated  from 
the  Sm02  determination.  Since  NIRS  cannot  distinguish  myoglobin 
absorbance  from  hemoglobin  absorbance,  Sm02  refers  to  combined 
saturation  of  hemoglobin  and  myoglobin  and  Pm02  refers  to  partial 
pressure  of  oxygen  in  the  interstitial  fluid.  The  same  absorption 
spectra  are  also  processed  by  the  computer  to  calculate  pH.  Subject- 
to-subject  differences  in  muscle  fiber  structure  produce  optical  vari¬ 
ations  that  are  corrected  prior  to  pH  calculation  (35).  The  corrected 
spectra  are  then  statistically  related  to  previously  determined  invasive 
measurements  of  interstitial  fluid  pH  to  produce  the  calibration  equa¬ 
tion  (28,  37).  This  second  calibration  equation  is  also  contained  in  the 
computer  and  used  to  calculate  pHm  from  the  corrected  spectra. 

The  light  output  of  the  system  is  assessed  with  three  NIST- 
traceable  reflectance  standards  (Avian  Technologies,  Wilmington, 
OH)  with  nominal  values  of  2,  50,  and  99%  prior  to  use  on  each 
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subject  to  allow  for  determination  of  the  absolute  values  of  SmC>2, 
PmCh,  and  pHm  (29). 

Statistical  analysis.  Values  for  NIRS  measurements  and  hemody¬ 
namic  parameters  are  presented  as  the  mean  ±  SE.  Data  from  all  the 
noninvasive  sensors  were  collected  continuously.  For  each  LBNP 
level,  the  last  3  min  of  data  for  each  measured  parameter  were 
averaged  to  provide  a  single  value  for  that  level. 

To  evaluate  the  overall  relationship  between  SV  and  TPR  and 
SmCh,  the  percentage  change  from  baseline  was  calculated  to  com¬ 
pensate  for  different  baseline  values  between  subjects.  Each  measured 
parameter  was  analyzed  using  a  linear  mixed  model  analysis  of 
variance  with  a  first  order  autoregressive  covariance  structure  to 
determine  if  there  was  a  significant  variation  during  progressive 
LBNP.  This  type  of  analysis  takes  into  account  the  repeated  nature  of 
the  experimental  design.  If  statistical  differences  were  found.  Bonfer- 
roni-corrected  comparisons  with  baseline  measurements  were  per¬ 
formed  to  determine  the  first  level  of  LBNP  that  could  be  distin¬ 
guished  statistically  from  baseline  (P  <  0.05).  Statistical  analysis  was 
performed  using  SPSS  (version  14.0,  SPSS,  Chicago,  IL). 

RESULTS 

All  10  subjects  completed  the  study.  Six  subjects  reached 
cardiovascular  collapse  at  —80  mmHg  LBNP  and  four  at  —90 
mmHg  LBNP.  In  all  subjects,  HR  returned  to  baseline  within 
1  min  after  the  LBNP  was  terminated.  Baseline  values  for 
specific  hemodynamic  and  NIRS-determined  metabolic  param¬ 
eters,  the  value  at  which  each  parameter  was  statistically 
different  from  baseline,  and  the  level  of  LBNP  at  which  this 
change  was  detected  are  presented  in  Table  1.  Arterial  oxygen 
saturation  was  not  altered  from  its  average  baseline  value  of 
98%  throughout  the  entire  LBNP  protocol.  MAP  was  only 
sensitive  to  the  detection  of  impending  cardiovascular  collapse 
in  that  it  essentially  signaled  the  end  of  the  LBNP  protocol.  HR 
changed  significantly  during  the  —60  mmHg  level,  while  PP 
decreased  significantly  during  the  —45  mmHg  level.  SV  was 
the  most  sensitive  of  the  hemodynamic  parameters,  with  a 
significant  decrease  occurring  at  the  first  level  of  LBNP  (—15 
mmHg).  Early  (—15  mmHg  LBNP)  metabolic  indicators  of 
hemodynamic  instability  included  PmCL  and  SmCL  deter¬ 
mined  by  NIRS.  In  contrast,  pH  did  not  decrease  from  baseline 
until  —70  mmHg. 

Figure  1  shows  the  average  (n  =  10)  of  each  of  the  measured 
and  calculated  parameters  for  every  level  of  LBNP  through 
—  80  mmHg.  HR  and  PP  were  unresponsive  until  —45  mmHg, 
when  HR  began  to  climb  and  PP  decreased.  SV,  PmCh,  and 
Sm02  responses  were  more  sensitive,  decreasing  from  the 
onset  of  LBNP  (—15  mmHg)  and  continuing  to  decrease  with 
the  application  of  greater  negative  pressure.  Muscle  pH  was 
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unaltered  until  —70  mmHg,  when  it  began  to  significantly  fall. 
The  overall  change  in  the  absolute  value  of  pHm  was  very 
small  (0.03  pH  units). 

Figure  2  illustrates  the  mean  percentage  change  from  base¬ 
line  for  SmCL  compared  with  SV  (Fig.  2A)  and  TPR  (Fig.  2 B) 
for  each  level  of  LBNP.  As  LBNP  progressed,  SmCL  declined 
in  parallel  with  SV  (r2  =  0.96)  and  inversely  with  TPR  (r2  = 
0.94).  PmCL  was  also  highly  correlated  with  SV  (r2  =  0.94) 
and  TPR  (r2  =  0.93). 

DISCUSSION 

The  major  finding  of  this  study  is  that  muscle  oxygen  (either 
SmCL  or  PmCL)  determined  noninvasively  with  a  novel  NIRS 
methodology  was  correlated  with  SV  and  was  the  earliest 
indicator  of  progressive  central  hypovolemia  in  human  volun¬ 
teers.  The  reductions  in  SV,  SmCL,  and  PmCL  were  apparent  at 
an  LBNP  level  that  represented  an  equivalent  blood  loss  of 
only  400-500  ml  (9).  Perhaps  more  importantly,  reductions  in 
muscle  oxygen  were  detected  much  earlier  than  the  standard 
clinical  measures  of  HR  and  BP.  A  further  finding  from  this 
study  is  that  a  small  but  significant  decrease  in  noninvasively 
measured  pHm  was  observed  near  the  same  time  as  an  increase 
in  HR,  but  later  than  a  detectable  change  in  PP. 

The  observation  that  decreases  in  SmCL  proportionately 
tracked  the  progressive  fall  in  SV  throughout  the  course  of  the 
LBNP  protocol  suggests  that  this  noninvasive  measurement  is 
a  sensitive  marker  of  a  limited  Q  and  delivery  of  blood  to 
peripheral  tissues  during  central  hypovolemia.  Likewise,  the 
inverse  linear  relationship  between  SmCL  and  TPR  implies  that 
skeletal  muscle  vasoconstriction  and  a  subsequent  reduction  in 
local  tissue  blood  flow  in  response  to  central  hypovolemia  is  a 
major  cause  of  reduced  regional  oxygen  in  this  laboratory 
model.  These  observations  are  consistent  with  the  work  of 
Fadel  et  al.  (1 1)  who  simultaneously  measured  relative  changes 
in  forearm  oxygenated  hemoglobin  plus  myoglobin,  blood 
flow,  and  vascular  conductance  during  LBNP  between  —10 
and  —50  mmHg.  These  authors  showed  strong  correlations 
between  NIRS  changes  and  both  blood  flow  and  vascular 
conductance.  Our  work  demonstrates  that  there  is  an  inverse 
relationship  between  NIRS-determined  tissue  oxygenation  and 
TPR  from  the  earliest  levels  of  hypovolemia  to  the  point  of 
presyncope. 

It  has  been  known  for  many  years  that  one  of  the  earliest 
compensatory  mechanisms  in  hemorrhage  is  sympathetically 
mediated  reflex  vasoconstriction  (1).  During  hemorrhage,  va¬ 
soconstriction  is  heterogeneous  throughout  the  body,  predom- 


Table  1.  Hemodynamic  and  metabolic  parameters 


Parameter 

Baseline  Value 

Statistically  Significant  Value 

LBNP  Level  (mmHg)  Where 
Statistically  Different  From  Baseline 

Arterial  oxygen  saturation,  % 

97.7±0.4 

NA 

none 

Mean  arterial  pressure,  mmHg 

101 .3  ±  3. 1 

9 1 .2  ±  5.3 

-90 

Heart  rate,  beats/min 

64.6±6.2 

86.8±5.6 

-60 

Pulse  pressure,  mmHg 

56.4±3.0 

49.2±2.5 

-45 

Stroke  volume,  ml 

140.2±7.4 

129.8±7.9 

-15 

NIRS  muscle  SO2,  % 

63.0±2.7 

58.2±2.7 

-15 

NIRS  muscle  PO2,  mmHg 

33.3  ±1.5 

30.8  ±  1 .3 

-15 

NIRS  muscle  pH,  pH  units 

7.44±0.01 

7.41  ±0.01 

-70 

Values  are  means  ±  SE  for  baseline  and  lower  body  negative  pressure  (LBNP)  level  of  significant  change.  NA,  not  applicable;  NIRS,  near  infrared 
spectroscopy. 
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Fig.  1.  Hemodynamic  and  near  infrared  spectroscopic  (NIRS)  derived  parameters  during  lower  body  negative  pressure  (LBNP).  Mean  ±  SE  for  all  subjects  (n  = 
10).  Negative  pressure  is  increased  progressively  from  baseline  (0).  *P  <  0.05  compared  with  baseline. 


inantly  affecting  skeletal  muscle  and  splanchnic  circulations  to 
redirect  blood  flow  to  the  heart  and  brain  (5,  13,  22).  Our 
results  complement  a  previous  suggestion  that  progressive 
LBNP  simulates  the  acute  hemodynamic  responses  occurring 
during  the  early  stages  of  hemorrhagic  shock  (9).  One  of  the 
novel  aspects  of  the  current  study  is  the  demonstration  that 
noninvasively  determined  skin  and  fat  corrected  Sm02  can  be 
used  as  an  indication  of  this  diversion  of  blood  from  the 
skeletal  muscle  (1),  suggesting  that  it  may  be  an  early  indica¬ 
tion  of  hemorrhagic  shock. 

Figure  1  shows  that  pHm  declined  slowly  until  there  was  a 
large  change  in  oxygen  delivery,  as  evidenced  by  considerable 
reductions  in  SV  and  Pm02.  pHm  is  significantly  decreased  at 
an  LBNP  level  of  —70  mmHg,  corresponding  to  a  fall  in  PmCL 
to  27.4  ±1.7  mmHg.  We  previously  demonstrated  a  similar 
relationship  between  tissue  pH  and  tissue  PCL  in  swine  liver 


during  hemorrhagic  shock  (27).  In  this  animal  study,  measure¬ 
ments  of  tissue  pH  and  PO2  with  invasive  sensors  demon¬ 
strated  a  breakpoint  (PCLcrit)  where  pH  began  to  decrease  after 
tissue  PO2  had  decreased  to  22.3  ±  3.8  mmHg  (linear  model; 
Ref.  27),  a  value  similar  to  that  observed  in  this  human  study. 
These  data  provide  evidence  that  LBNP  causes  a  tissue  meta¬ 
bolic  disturbance  (i.e.,  muscle  dysoxia  in  the  forearm)  similar 
to  that  observed  during  actual  bleeding. 

PP  has  been  suggested  as  an  easy  to  obtain  hemodynamic 
parameter  that  provides  an  early  indication  of  the  reduction  in 
SV  and  oxygen  delivery  (6).  In  a  study  of  30  trauma  patients 
requiring  helicopter  transport,  PP  was  a  better  predictor  of 
mortality  than  MAP,  HR,  or  arterial  oxygen  saturation  (8).  The 
current  study  shows  that  noninvasive  Pm02  is  an  even  earlier 
indicator  of  reduced  oxygen  delivery  than  PP  (Table  1,  Fig.  1). 
Likewise,  transcutaneous  oxygen  tension  (PtcC>2)  has  been 


A 

aT 

D) 

C 

03 

Fig.  2.  Percent  change  from  baseline  for  o 
each  level  of  LBNP.  NIRS  skin  and  fat  ^ 
corrected  deep  muscle  oxygen  saturation  ^ 
(SmC>2)  is  evaluated  as  a  function  of  stroke  Q 
volume,  r2  =  0.96  (A)  and  total  peripheral  E 
resistance,  r2  =  0.94  ( B ).  Error  bars  repre-  ^ 
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assessed  in  emergency  patients  (31);  nonsurvivors  had  signif¬ 
icantly  lower  PtcC>2  values  compared  with  survivors.  However, 
application  of  PtcCL  equipment  to  rapidly  deteriorating  pa¬ 
tients  is  limited  by  the  necessary  20  min  equilibration  time 
needed  to  heat  the  skin  (31).  Such  a  heating  protocol  is  not 
necessary  for  NIRS  technology,  suggesting  that  it  may  be  more 
appropriate  for  monitoring  trauma  patients.  Indeed,  the  spec¬ 
troscopic  assessment  of  SmCL  and  Pm02  is  noninvasive,  con¬ 
tinuous,  portable,  and  can  be  started  rapidly  in  a  trauma  or 
emergency  response.  This  monitoring  technique  has  the  poten¬ 
tial  to  provide  early  warning  of  impending  cardiovascular 
collapse  in  the  pre-  and  early  hospital  settings. 

Our  laboratory  model  of  progressive,  central  hypovolemia 
(LBNP)  provides  reproducible  hemodynamic  responses  both 
within  and  between  subjects;  our  hemodynamic  results  were 
similar  to  those  obtained  in  previous  studies  in  the  same  and 
other  laboratories  (7,  9,  23).  While  this  model  provides  us  with 
some  insight  into  changes  in  tissue  oxygenation  during  the 
initial  phase  of  hemorrhagic  shock,  we  cannot  extrapolate  our 
findings  to  monitoring  of  patients  in  septic  shock.  The  micro- 
vascular  response  to  sepsis  is  expected  to  have  a  different  time 
course  and  oxygenation  profile  and  is  the  subject  of  a  separate 
patient  study. 

The  NIRS  techniques  used  in  this  study  are  unique  in  that 
they  report  values  for  PmCL,  SmCL,  and  pHm  deep  within 
the  muscle,  corrected  for  light  absorption  by  skin  and  fat. 
Other  studies  employing  NIRS  report  only  tissue  oxygen 
saturation,  where  the  measurement  is  a  nonspecific  weighted 
average  of  oxygen  saturation  in  the  skin,  fat,  and  muscle 
depending  on  the  design  of  the  sensor  (3,  10).  Our  sensor  is 
specifically  designed  to  collect  and  analyze  spectra  from  the 
muscle,  eliminating  spectral  interferences  from  skin  pig¬ 
ment  and  fat  (34).  In  addition,  the  calculation  methods  used 
to  determine  Pm02  and  Sm02  are  independent  of  variations 
in  muscle  optical  properties  between  subjects.  However,  the 
pHm  calculation  methods  rely  on  applying  a  correction 
factor  to  account  for  subject-to-subject  variation  in  muscle 
optical  properties  (35).  The  robustness  of  this  correction 
factor  may  limit  the  absolute  accuracy  of  the  pHm  measure¬ 
ment  for  sick  individuals.  Current  work  is  investigating  the 
limitations  of  the  pH  calculation  and  in  providing  accurate 
pHm  determinations  in  both  healthy  and  sick  individuals. 
The  Pm02  calculation  described  in  the  appendix  does  not 
correct  the  hemoglobin  dissociation  curve  for  changes  in 
pH,  however,  for  in  this  study,  pH  changes  are  very  small. 
Future  work  will  enhance  the  PmC^  assessment  by  incor¬ 
porating  this  additional  information. 

Conclusions.  We  demonstrated  a  noninvasive  technique 
for  the  determination  of  Sm02,  PmCL,  and  pHm  in  a  clinical 
laboratory  model  that  simulates  hemorrhage  in  humans. 
Sm02  and  PmCH  were  found  to  be  early  indicators  of  a 
reduction  in  oxygen  delivery,  most  likely  resulting  from 
increased  vasoconstriction  to  shunt  blood  from  the  skeletal 
muscle  to  the  heart  and  brain.  Muscle  pH  begins  to  decrease 
only  when  oxygen  levels  are  significantly  reduced,  suggest¬ 
ing  that  the  decrease  in  pHm  may  signal  the  onset  of  dysoxia 
in  the  muscle  tissue.  This  noninvasive,  continuous  technique 
lends  itself  to  monitoring  critically  ill  patients  at  risk  for 
hemodynamic  instability. 


APPENDIX 


The  method  for  determining  muscle  oxygen  saturation  (S1T1O2)  is 
reported  in  detail  by  Yang  et  al.  (36),  but  is  described  briefly  here. 
SmCH  is  defined  by  Eq.  1: 

C(Hb02+Mb02)<T> 

S  m  (_)  2  — 

C(Hb02+Mb02)<'T>  +  C(Hb+Mb)<T>  (7) 

_  QtibOj+MKy 

C(Hb02+Mb02)  +  C(Hb+Mb) 


where  L  is  the  average  photon  pathlength  through  the  tissue, 
k-Hbo2+Mbo2  is  the  oxygenated  heme  concentration,  and  CHb+Mb  is  the 
deoxygenated  heme  concentration,  Hb  represents  hemoglobin  and  Mb 
myoglobin. 

Attenuation  of  light  by  the  subject’s  tissue,  Aexp(X)  at  wavelength 
X  is  defined  in  Eq.  2: 


Aexp(k) 


In 


Iref(k) 

.  HX) 


(2) 


where  Iref(k)  is  the  measured  diffuse  reflectance  intensity  from  a  99% 
diffuse  reflectance  reference  standard  at  wavelength  X,  I(X)  is  the 
measured  diffuse  reflectance  intensity  from  the  subject  at  wavelength 
X.  Incident  light  attenuation  by  tissue  is  caused  by  both  absorption 
and  scattering  events.  Light  is  absorbed  by  hemoglobin  in  the  small 
blood  vessels  and  myoglobin  in  the  cells,  as  well  as  both  intravas¬ 
cular  and  extravascular  water  and  melanin  pigment  in  the  skin. 
Light  is  scattered  away  from  physical  structures  in  the  tissue  such  as 
blood  vessels  and  muscle  fibers,  as  well  as  fat  that  lies  over  the 
muscle. 

Muscle  or  tissue  oxygen  saturation  can  be  calculated  from  mea¬ 
sured  spectra  by  assuming  that  light  attenuation  is  primarily  a  result  of 
absorption  by  oxygenated  and  deoxygenated  hemoglobin  and  myo¬ 
globin,  as  well  as  a  result  of  the  scattering  by  the  tissue.  Several 
methods  have  been  proposed  for  acquiring  the  spectra  and  performing 
the  calculations;  these  methods  are  reviewed  for  application  in  phys¬ 
iology  by  Ferrari  et  al.  (12).  All  of  these  methods  involve  fitting  the 
measured  attenuation  spectrum  to  different  mathematical  models  that 
described  light  penetration  through  the  tissue  to  separate  absorption 
from  scattering  and  then  fitting  the  obtained  absorption  spectra  to  the 
Beer-Lambert  Law,  which  relates  the  absorption  to  the  heme  concen¬ 
tration.  Primary  differences  between  methods  relate  to  how  light 
scattering  is  measured  and  separated  from  the  absorption.  Light 
scattering  can  be  measured  directly  with  spatially  resolved  continuous 


Oxygenation  Index 


a> 


LBNP  (mm  Hg) 

Fig.  3.  Change  from  baseline  oxygenation  index  (01).  OI  at  baseline  is  set  to 
zero,  changes  measured  in  optical  density  (OD)  units.  Mean  ±  SE  for  all 
subjects  («  =  10).  Negative  pressure  is  increased  progressively  from  baseline 
(0).  *P  <  0.05  compared  with  baseline. 
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wave  (CW)  spectroscopy,  time-resolved  spectroscopy,  and  frequency- 
domain  spectroscopic  techniques.  Recently  van  Beekvelt  et  al.  (32) 
and  Homma  et  al.  (14)  demonstrated  that  fat  layers  can  influence  the 
accuracy  of  an  SmC>2  calculation.  We  demonstrated  a  method  to 
correct  attenuation  spectra  for  the  spectral  interference  of  fat  prior  to 
their  use  in  calculating  physiologically  important  parameters  (34). 

To  calculate  Sm02,  we  first  remove  the  components  of  the  spec¬ 
trum  that  result  from  skin  pigment  absorption  and  fat  scattering  by 
using  a  two-source  fiber  optic  probe.  One  source  is  placed  close  (2.5 
mm)  to  the  fiber  optic  bundle  that  transmits  light  to  the  spectrometer. 
This  captures  light  from  only  the  skin  and  fat  layers.  The  second 
source,  a  farther  distance  from  the  detector  bundle  (30  mm),  captures 
light  from  the  skin,  fat,  and  muscle  layers.  Light  collected  from  the 
short  distance  pair  is  orthogonalized  with  the  light  from  the  long 
distance  pair  to  generate  a  spectrum  that  describes  attenuation  from 
only  the  muscle  layer.  The  details  of  this  method  are  described  in  a 
prior  publication  (34).  The  corrected  spectrum  is  then  used  for  the 
calculation  of  SmCL  based  on  an  adaptation  of  a  method  first  proposed 
by  Stratonnikov  and  Loschenov  (30).  In  this  method,  a  Taylor  expan¬ 
sion  attenuation  model  is  used  to  model  light  that  penetrates  through 
the  tissue.  The  light  absorption  is  modeled  by  Beer’s  law  and  that  of 
the  scattering  is  modeled  as  a  Taylor  expansion  term. 

The  skin  color  and  fat  corrected  attenuation  spectrum  is  then 
described  with  the  Taylor  series  expansion  model  Am0dei  detailed  by 
Stratonnikov  and  Loschenov  (30) 

/I0(X)\ 

Amodel(^-)  —  ln  _  (C°  Cl^  f  CHb+Mb  e(^-)  .j. 

+  cHb02+Mb02£Hb02  (V)  +  ^wat  ^wat  (X)}  In  (10) 

where  Io  is  the  light  source  intensity,  Co  and  Ci  are  constants,  cwat  is 
the  concentration  of  water.  €hi>,  eHbc>2  •  and  ewat  are  known  extinction 
coefficients  of  Hb,  HbCL,  and  water,  respectively.  Since  hemoglobin 
and  myoglobin  have  nearly  identical  extinction  coefficients,  only  the 
extinction  coefficients  of  hemoglobin  are  required.  The  function  (co  + 
Ci\)  describes  the  portion  of  the  spectrum  resulting  from  light  that  is 
scattered,  as  well  as  the  wavelength-independent  absorption  caused  by 
chromophores  other  than  heme  and  water  and  the  experimental  dif¬ 
ference  caused  by  using  Iref  in  Eq.  2  as  the  incident  light  intensity  (Io) 
when  calculating  the  attenuation  spectra. 

c0,  ci  as  well  as  cHb+Mb,  cHbo2+Mbo2  ,  cwat,  and  L  are  obtained  by 
nonlinear  least  square  fitting  of  the  measured  attenuation  spectrum 
Aexp  to  the  modeled  spectrum  Am0dei  described  by  Eq.  3.  Once 
CHb+Mb  and  cHb02+Mbo,  are  obtained,  Sm02  is  calculated  using  Eq.  1. 
One  of  the  limitations  of  this  technique  is  the  inability  to  separate  the 
deoxygenation  of  myoglobin  from  hemoglobin,  hence,  the  general 
name  muscle  oxygen  saturation. 

The  partial  pressure  of  oxygen  in  the  muscle  (Pm02)  can  be 
calculated  from  Sm02  determined  with  this  method.  Severinghaus 
(25)  describes  the  relationship  between  SO2  and  PO2  under  standard 
physiological  conditions: 

P02  =  exp(0.385-ln  (S02~‘  -  l)-1  +  3.32  (4) 

-  (72-S02r‘  -  (S02)/6) 

For  our  application,  Sm02  is  calculated  from  skin  pigment  and 
fat-corrected  spectra.  Then  Pm02  is  calculated  using  Eq.  4.  Again, 
Pm02  represents  an  average  of  hemoglobin  and  myoglobin  oxygen 
tension.  We  also  recognize  that  calculation  of  Pm02  should  be 
compensated  for  variation  in  pH  and  PC02.  Figure  1  shows  that 
changes  in  pH  were  minimal  in  this  LBNP  model,  with  significant 
decreases  in  pH  occurring  late  in  the  progression.  For  this  reason, 
measurement  of  Pm02  during  early  hypovolemia  will  be  unaffected 
by  this  omission. 

The  calculations  were  implemented  in  software  written  in  Matlab 
[version  7.0.4.365  (R14)  Service  Pack  2,  Mathworks,  Natick,  MA], 
The  fitting  algorithm  was  implemented  with  the  “lsqcurvefit”  function 
in  the  Matlab  optimization  toolbox  version  3.0.2. 


Comparison  with  alternate  CW  NIRS  techniques.  The  new  NIRS 
methodology  is  compared  with  a  traditional  analysis  technique.  Spec¬ 
tra  from  the  long  distance  sensor,  prior  to  correction  for  skin  color  and 
fat,  were  used  with  the  UCL6  algorithm  (19).  This  algorithm  has  been 
used  on  the  NIRO-1000  instrument  (Hamamatsu  Photonics  KK).  The 
algorithm  calculates  relative  values  of  oxyheme  (AHb02)  and  deoxy- 
heme  (AHb)  with  six  discrete  values  from  the  NIR  spectrum.  Delta 
represents  the  change  from  baseline  absorbance  at  the  six  key  wave¬ 
lengths  in  optical  density  units.  To  compare  to  our  calculated  Sm02, 
we  then  used  the  AHb02  and  AHb  values  to  determine  the  oxygen¬ 
ation  index  (OI)  as  described  in  Ferrari  et  al.  (12). 

OI  =  AHb02  -  AHb  (5) 

OI,  a  measure  of  the  relative  changes  in  tissue  oxygenation  is 
illustrated  in  Fig.  3.  This  relative  measure  shows  the  same  trends  and 
statistical  significance  as  the  determination  of  Pm02  (Fig.  ID)  and 
Sm02  (Fig.  1£).  OI  does  not  take  into  account  subject  variation  in 
skin  color,  fat  content,  or  muscle  scattering,  so  only  relative  changes 
in  tissue  oxygenation  during  progressive  LBNP  can  be  assessed  (19). 

Our  study  uses  LBNP  to  simulate  progressive  bleeding  and  seeks  to 
identify  a  noninvasive  spectroscopic  parameter  that  can  be  used  by 
first  responders  to  determine  if  a  patient  is  at  risk  for  cardiovascular 
collapse.  Since  it  will  often  be  the  case  that  the  baseline  condition  will 
be  unknown,  relative  measures  of  tissue  oxygenation  will  be  of 
limited  value.  The  new  NIRS  algorithms  presented  here  produce 
repeatable  values  for  Sm02  and  Pm02  across  multiple  subjects  at 
varying  levels  of  hypovolemia-induced  vasoconstriction.  These  pa¬ 
rameters  have  the  potential  to  be  useful  in  the  treatment  of  critically 
ill  patients;  the  next  step  is  to  validate  the  technique  in  actual  patients. 
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